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Abstract 

The hadronic decays of the r lepton can be used to determine the effective 
charge ar{''n'^i) for a hypothetical r-lepton with mass in the range < m^/ < 
rrir. This definition provides a fundamental definition of the QCD coupling 
at low mass scales. We study the behavior of cxr at low mass scales directly 
from first principles and without any renormalization-scheme dependence by 
looking at the experimental data from the OPAL Collaboration. The results 
are consistent with the freezing of the physical coupling at mass scales s = rn^, 
of order 1 GeV^ with a magnitude a,- ~ 0.9 it 0.1. 
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I. INTRODUCTION 



One of the major uncertainties in making reliable predictions in quantum chromodynam- 
ics (QCD) is to understand the theory at low momentum scales where the coupling becomes 
large and non-perturbative effects become important. In fact, it is well known that pertur- 
bation theory itself is not well defined in the infrared domain, since the perturbative series is 
asymptotic. For infrared-safe quantities the non-perturbative effects can be parameterized 
as power corrections of the form c„(A/(5)", where A is the QCD-scale and Q is the hard 
scale of the process considered. The coefficients c„ cannot be calculated in general, but in 
some cases such as the event-shapes in e'^e" annihilation, one can determine the dominating 
power n and also find relations between the c„ for different observables by using the operator 
product expansion [1] or renormalon calculus (for a review see ref. [2]). 

The behavior of fixed-order perturbation theory at low momentum scales is governed by 
the running coupling asik!^)- The conventional coupling q;ms(A;^), defined using the modified 
minimal subtraction renormalization scheme and dimensional regularization, is analytically 
singular at a scale A;^ = A|jg. This nonphysical behavior leads to a number of difficulties, 
including renormalon n\ growth of the coefficients when one makes perturbative expansions 
of physical observables in the MS scheme. These problems can be traced to the fact that 
integrals over the running coupling which appear in bubble graphs are ill-defined due to the 
non-analyticity of the MS coupling. 

There also exists the possibility to define a running coupling which stays finite in the in- 
frared. One such example is the "time-like" effective coupling which is used in the dispersive 
approach [3-5]. In this approach, an observable is written as an integral over a process de- 
pendent characteristic function times a universal "time-like" effective coupling. The idea is 
that such a coupling may give an effective measure of the interaction at low scales [5] . Power 
corrections can then be written as process dependent moments over the effective coupling. 

An alternative procedure is to define the fundamental coupling of QCD from a given 
physical observable [6,7]. These couplings, called effective charges, are all-order resumma- 
tions of perturbation theory and include all non-perturbative effects. Since these physical 
charges correspond to the complete theory of QCD, it is guaranteed that they are analytic 
and non-singular. For example, it has been shown that unlike the MS coupling, a physical 
coupling is analytic across quark flavor thresholds [8,9]. Furthermore, we expect that a 
physical coupling should stay flnite in the infrared when the momentum scale goes to zero. 
In turn, this means that integrals over the running coupling are well defined for physical 
couplings. An additional question is whether the physical couplings freeze to a constant 
value in the infrared. 

Once such a physical coupling a^hysik"^) is chosen, other physical quantities can be ex- 
pressed as expansions in a^nys by eliminating the MS coupling which now becomes only an 
intermediary [10]. In such a procedure there are in principle no further renormalization 
scale (//) or scheme ambiguities. The physical couplings satisfy the standard renormaliza- 
tion group equation for its logarithmic derivative, dccphys/dln A;^ = /3phys[cKphys(^^)]) where 
the first two terms in the perturbative expansion of the Cell-Mann Low function /Sphys are 
scheme-independent at leading twist whereas the higher order terms have to be calculated 
for each observable separately using perturbation theory. 



2 



Quantum field theoretic predictions which relate physical observables cannot depend on 
theoretical conventions such as the choice of renormalization scheme or scale (/i). The most 
well-known example is the perturbative "generalized Crewther relation" [11] in which the 
leading twist QCD corrections to the Bjorken sum rule for polarized deep inelastic scattering 
at a given lepton momentum transfer are related through a geometric series to the QCD 
corrections to Re+e- a corresponding CM energy squared, s* = s*{Q'^), independent of 
renormalization scheme, [l+aR{s*) /7i){l—ag^{Q'^) /n) = 1 [12]. The ratio of the scales 
has been computed to NLO in PQCD. Such leading- twist predictions between observables 
are called "commensurate scale relations" and are identical for conformal and nonconformal 
theories [10]. In addition, the conformal coefficients are free of the renormalon factorial 
growth [13-15] 

For example, in QED, the Cell-Mann Low running coupling q;qed(^^) = i-u{k^k^) ' "^^ich 
is formally defined from the renormalization of the dressed photon propagator, is a physical 
coupling since it could be determined from a measurement of the part of the potential V{k'^) 
between two infinitely heavy test charges which is linear in their charges; i.e., 47rQ;QED(^^) = 
— /c^Viin(fc^). Using the skeleton expansion [16], the coefficients in perturbative expansions 
of other physical quantities are identical to that in a theory which is conformal, since all 
effects of the non-zero /J function are already summed into the integrals over the running 
coupling. By the mean value theorem, the same is also true for the standard perturbative 
expansion if the scale, at which to evaluate the coupling, is properly chosen [17,18]. 

It is not as simple to identify a suitable physical coupling to be used in the case of 
QCD. For a skeleton expansion to be possible, the Abelian part of the coupling should 
coincide with the Gell-Mann Low coupling, since QCD becomes an Abelian theory in the 
analytic limit A^^c ^ at fixed Cpc^s and fixed Uf/Cp [19]. A possible candidate for a 
physical coupling in QCD which fulfills this requirement is the av{k'^) scheme defined from 
the potential between two heavy test color charges [20]. However, in contrast to the Abehan 
case, this definition is problematic since the "if— graphs" which arise from gluon exchange 
diagrams with a horizontal gluon rung connecting the "first" and "last" exchanged gluons 
have an infrared sensitivity which depends on the details of the test charge wavefunction [21]. 
Another possible generalization of the Cell-Man Low coupling to non-Abelian theories is 
the "pinch" scheme [22-24] which rearranges the contributions to scattering amplitudes to 
insure a structure Q;pinch(A;^) — i-nlk^'^k^^) similar to that of QED. As in QED, expansions 
in the pinch scheme have the same structure as those of a conformal theory. The pinch 
charge is a promising physical scheme, but at this time the complexity of higher order 
calculations in this scheme and its indirect connection to measurements has prevented its 
practical implementation, although recently there has been attempts to make an all-orders 
definition of a QCD effective charge [25] . 

In this note we will discuss an alternative definition of a physical coupling for QCD 
which has a direct relation to high precision measurements of the hadronic decay channels 
of the T~ — >■ i/^h~. Details on the extraction of as^ni^) from r decays can be found in 
refs. [4,26-35] (for some recent developments in the extraction of as{ml) from r decays see 
for example [36-39]). 

^et Rt be the ratio of the hadronic decay rate to the leptonic one. Then R^ = 



Rl 



1 + 



where i?^ is the zeroth order QCD prediction, defines the effective charge 
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ttr- Throughout this paper wc will concentrate on non-strange decay modes and thus 
= 3S'Ew|KidP where S'ew = 1-0194 is an electroweak correction term [40] and |KidP = 
0.9512 ± 0.0008 [41] is the relevant CKM matrix element. The data for r decays is well- 
understood channel by channel, thus allowing a precise separation of vector and axial- vector 
decay modes which can therefore be studied separately. 

The measured invariant mass spectrum for the non-strange hadronic decay modes can 
also be used to study hypothetical r-leptons with a smaller mass, m^/ < rrir [33-35]. In this 
way the r-decay data allows us to study the behavior of the coupling q;^(s) in the region 
< s < ml and address the question whether this physical coupling freezes in the infrared. 



II. EXTRACTING a^-iS) FROM DATA ON r-DECAYS 

The experimental data on the non-strange hadronic r-decays can be used to define the 
hadronic decay rate normalized to the leptonic one for a hypothetical r-lepton with mass in 
the range < m,-' < mr in the following way [34]: 



m^ 



127r5'Ew|Kid| 



ds 



mz 



2 r 



mz 



(0) 



r(i) 



where 27ilmllyjj^(s] 



v/a{s) are the spectral functions for the measured non-strange 

hadronic final states with angular momentum J = 1. The scalar contribution Imny^^(s) is 
assumed to vanish for the vector current and is given by the single pion pole for the axial 
current. The above definition coincides with the standard definition in the case m-r' — m^- 
Note, however, that the right-hand side above uses m-r' instead of m-r not only in the upper 
integration limit but also in the kinematic prefactors. This way the end-point is suppressed 
in the same way for the hypothetical r-lepton as for the real one [34]. The ratio R^'^{ml,) 
is thus defined as if a hypothetical lepton of mass m,-' existed. 

Since states with non-zero strangeness can be excluded, all of the hadrons in the final 
state can be assumed to arise from t~ — > VtVl<\. We can then define effective charges 
oi^/^{ml,) for the vector and axial- vector decay modes as follows: 
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1 + 



a 



\'/A 



TT 



(2) 



Notice that the combination — does not receive a perturbative QCD contribution 
at leading twist; thus if QCD is correct this contribution should be power law suppressed 
at high energies. It is thus natural to identify the complimentary R^ + R^ combination 
which has canonical perturbative QCD contributions as the preferred QCD effective charge, 
aT-{ml,), defined by 



Rriml,) = R^{ml,) + R^im^,) = R° 



1 + 



a-rimz 



TT 



(3) 



For completeness, we also recall how one can use experimental data to measure the decay 
ratio of hypothetical r-leptons for masses well above m^- [42]. Just as in the case of r-leptons 
one can define a local unintegrated effective charge aR{s) directly from the annihilation data: 
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i?e+e- (s) 



(j(e+( 
(7(e^ 



hadrons) 
//+//-) 



1 + 



TT 



(4) 



where -Rg+g- is the zeroth order QCD prediction. If we assume isospin invariance, then the 
decay ratio of a hypothetical r-lepton in the vector channel R^{m^,) can be written as a 
spectral integral with weight f{x) = (1 — + 2x), where x = s/m^,, of the annihilation 

cross section e+e~ ^ 7* — ^ hadrons in the isospin 1 = 1 channel. This allows the measure- 
ment of aY (ml 



a 



(1=1) 

R 



well above the physical mass of the r [42] . We thus can relate the and 



effective charges: 



a 



V 



{ml.) 



Jo 



ds 



mt 



f 



mz 



Jo 



ds 



mz 



f 



a 



(1=1) I 

R y 



mz 



(5) 



The mean value theorem then implies 



a}:iml,)=at'\sn, 



< s* < mlr, 



(6) 



a form of commensurate scale relation. In the case of three flavors (nf = 3) the above relation 
is still valid to next-to-next-to-leading order even if one does not restrict oneself to the 1=1 
channel. To next-to- leading order in the scale s* is given by (see for example [10]) 



s = m^i exp 



19 169/^ 2 \aR 

11 Uf — 

12 576 V 3 7 TT . 



(7) 



Before continuing we also note that as an alternative definition of a hypothetical r-lepton 
with mass above m^- one could use ImnY''(s) measured from r-decays for the integration 
region < s < and data from i?e+e- (s) in the 1=1 channel for the remaining integration 
region m^ < s < m^.,. 

The empirical behavior of the decay ratios Rr, , R^, and R^ — R^ using data on r- 
decays as determined by the OPAL collaboration at LEP [35,43] are shown in fig. 1. There 
are several striking features: 

1. the i?y — R^ combination tends to vanish at high scales, showing that only higher 
dimensional operators contribute. This is a highly non-trivial test of QCD. [44,45] 

2. The Rr contribution has only a slow variation in the low mass range. 



are 



The corresponding effective charges a^, , and as well as the difference a 
shown in Fig. 2. 

Based on the analysis by the OPAL collaboration [35], the experimental value of the 
couphng ar{s) = 0.621 ± 0.008 at s = m^ corresponds to a value of aj^{M'^) = (0.117- 
0.122) ±0.002, where the range corresponds to three different perturbative methods used in 
analyzing the data. This result is, at least for the fixed order and renormalon resummation 
methods, in good agreement with the world average aMs(M|) = 0.117± 0.002 [46]. However, 
from the figure we also see that the effective charge only reaches ar{s) ~ 0.9 ± 0.1 at 
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s (GeV) 



FIG. 1. The non-strange hadronic decay rate of a hypothetical r lepton with m^, = s versus the 
upper integration limit s for the currents V , A, V + A, and V — A. Error bands include statistical 
and systematic errors. 



s = 1 GeV^, and it even stays within the same range down to s ~ 0.5 GeV^. This resuh is in 
good agreement with the estimate of Mattingly and Stevenson [47] for the effective couphng 
CiR{s) ~ 0.85 for ^/s < 0.3 GeV determined from e"'"e~ annihilation, especially if one takes 
into account the perturbative commensurate scale relation, a^(m^;) = ur^s*) where, for 
an = 0.85, we have s* ~ O.lOm^, according to Eq. (7). As we will show in more detail in 
the next section, this behavior is not consistent with the coupling having a Landau pole but 
rather shows that the physical coupling is much more constant at low scales, suggesting that 
physical QCD couplings are effectively constant or "frozen" at low scales. 

At the same time, it should be recognized that the behavior of ar{s) in the region 
s < 1 GeV^ is more and more influenced by non-perturbative effects as the scale is lowered. 
Even though the dominant non-perturbative effects cancel in the sum of the vector and 
axial-vector contributions as can be seen by looking at the corresponding effective charges 
individually. Looking at a^{s), we see that it more or less vanishes as the integration region 
moves to the left of the two-pion peak in the hadronic spectrum. In the same way the 
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FIG. 2. The effective charges for non-strange hadronic decays of a hypothetical r lepton with 
m^i = s versus the upper integration limit s for the currents V , A, V + A. Also shown is the 
difference of the effective charges for the vector and axial-vector current. Error hands include 
statistical and systematic errors. 



behavior of Oi^{s) at small scales is governed by the single pion pole. 



III. ANALYSIS OF THE INFRARED BEHAVIOR OF a^{S) 

In order to be able to analyze the infrared behavior of the effective coupling ar{s) in 
more detail, we will compare with (a) the fixed-order perturbative evolution of the ar{s) 
coupling on the one hand, and (b) with the evolution of couplings that have non-perturbative 
or all-order resummations included in their definition. For the latter case, many different 
schemes have been suggested, and we will concentrate on two of them: the one-loop "time- 
like" effective coupling aefr(s) [3-5], and the modified ay coupling calculated from the static 
quark potential using perturbative gluon condensate dynamics [48]. 

The perturbative couplings evolve according to the standard evolution equation 
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das(g) 
dins 



(8) 



where as(s) = as{s) / (Ati) . The first two terms in the /3-function, /3o and are universal at 
leading twist whereas the higher order terms are scheme dependent. Currently the /9-function 
is known to four loops {^s) in the MS scheme and to three loops (/32) in the ar scheme. In 
the latter case there also exists an estimate of the four-loop term. For completeness these 
terms are summarized in the appendix. 

Fig. 3 shows a comparison of the experimentally determined effective charge Q!r(s) with 
solutions to the evolution equation (8) for ar at two-, three-, and four-loop order normalized 
at m^. It is clear from the figure that the data on a^(s) docs not have the same behavior as 
the solution of the (universal) two- loop equation which is singular^ at the scale s ~ 1 GeV^. 
However, at three loops the behavior of the perturbative solution drastically changes, and 
instead of diverging, it freezes to a value a,- ~ 2 in the infrared. The reason for this 
fundamental change is, of course, the negative sign of f3r,2- At the same time, it must 
be kept in mind that this result is not perturbatively stable since the evolution of the 
coupling is governed by the highest order term. This is illustrated by the widely different 
results obtained for three different values of the unknown four loop term i3t-,3 which are also 
shown^ . Still, it may be more than a mere coincidence that the three-loop solution freezes in 
the infrared. Recently it has been argued that ^^(s) freezes perturbatively to all orders [49]. 
Given the commensurate scale relation (6) this should also be true perturbatively for aT-{s). 
It is also interesting to note that the central four-loop solution is in good agreement with 
the data all the way down to s ~ 1 GeV^. 

The one-loop "time-like" effective coupling [3-5] 



/ N 47r fl 1 

aes{s) = — <- arctan 



TT A2 



(9) 



is obtained from the analytic continuation of the one-loop coupling «s(Q^) = 
(47r)/[/3o ln((5^/A^)] which defines the spectral density dcteffl-s) /dins = [as{—s+ie)—as{—s — 
ie)]/{2TTi) in the same way as Re+e-{s) is related to the Adler D-function. The resulting 
effective coupling is finite in the infrared and freezes to the value aesis) — 47r//3o as s — > 0. 
It is also instructive to expand the "time-like" effective coupling for large s. 



An 



/9oln (s/A2) 



1 - 



1 



TT 



+ 



TT 



31n"(s/A2) 51n*(s/A2) 



+ ... 



3 I 47r 



The same divergent behavior would also be seen at three- and four- loop order in the MS scheme 
where both (32 and /?3 are positive for rif = 3. 

^The values of Pr,3 used are obtained from the estimate of the four loop term in the perturbative 
series of Rr, Kf^ = 25 ± 50 [30]. 
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FIG. 3. The effective charge ar for non-strange hadronic decays of a hypothetical t lepton with 
m^i = s compared to solutions of the fixed order evolution equation (8) for at two-, three-, and 
four-loop order. Error hands include statistical and systematic errors. 

This shows that the "time-hke" effective couphng is a resummation of {'k'^ l3Qa^)'^-corieci\OTa.s 
to the "space-hke" couphng which occurs in the analytic continuation. The "time-hke" 
effective couphng can also be defined to higher orders (for a recent review see [50]) but the 
infrared behavior persists also in these cases. 

The form for the finite coupling aefr given in Eq. (9) obeys standard PQCD evolution at 
LO. Thus one can have a solution for the perturbative running of the QCD coupling which 
obeys asymptotic freedom but does not have a Landau singularity. 

The evolution of the modified ay(Q^,Ag^) coupling is to leading order governed by the 
evolution equation, 



dav(Q^A^J 
dlng2 



-/?o(Agc/g)4(Q',A^c), 



where 



/5o(V/« = 5-12§+(l2|-6^ 



3 g \ Q + 2Age 2 
+ — I In — ;7T^ nf. 



2A 



g-2A 



(10) 



(11) 
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The difference compared to the ordinary evolution equation is thus that the gluonic contri- 
bution to Po freezes out at a scale Q ~ Age - the scale of the gluon condensate. 
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FIG. 4. The effective charge for non-strange hadronic decays of a hypothetical r lepton 
with m^, = s compared to the one-loop "time-like" effective coupling aeff(s) and the modified ay 
coupling calculated from the static quark potential using perturbative gluon condensate dynamics. 
Error bands include statistical and systematic errors. 

Fig. 4 shows a comparison of the experimentally determined effective charge ar{s) with 
the one-loop "time-like" effective coupling Q;efr(s) and the modified ay coupling calculated 
from the static quark potential using perturbative gluon condensate dynamics. In both 
cases the solutions have been normalized to the data at m,-. In addition, when solving 
the evolution equation for ay (10) we have used the value Age = 581 MeV which makes 
the solution to the one-loop evolution equation (8) agree with the data at rrir if one sets 
A = Age. As can be seen from the figure, the data on ar{s) agrees well qualitatively with 
both of these simple examples of freezing couplings down to the scale s ~ 1 GeV^. Below 
this scale the non-perturbative effects in ar{s) associated with the single pion pole in the 
axial vector current and the double pion peak in the vector current starts to dominate and 
makes a direct comparison with models which do not contain the spectrum of hadrons less 
meaningful. 
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IV. CONCLUSIONS 



The results of this paper show the advantages of defining the QCD couphng directly from 
a physical observable. The resulting effective coupling can be defined even at low scales, 
it is finite and analytic, and it has no scheme or renormalization scale ambiguities. As we 
have shown, the hadronic decays of the r lepton can be used to determine the effective 
charge «^(m^,) for a hypothetical r-lepton with mass in the range < m^/ < rrir. The 
high precision of r hadronic decay data thus provides a precise standard definition for the 
fundamental QCD couphng. QCD predictions for other observables can then be expressed 
as functions of this coupling 0(t-{s), thus relating observable to observable. 

An important feature of ar{s) is its apparent near constant behavior at low mass scales. 
The empirical results using the OPAL data are consistent with the freezing of the physical 
coupling at mass scales s = m^/ of order 1 GeV^ with a magnitude o;^ ~ 0.9 ± 0.1. These 
results contrast dramatically with standard expectations of a divergent coupling based on the 
universal two-loop coupling which becomes infinite at small mass scales. At higher orders of 
perturbation theory for the beta function, the behavior of the coupling is scheme dependent 
and in addition it is dominated by the highest order term. In the physical a^-scheme the 
behavior of the coupling in the infrared is not so clear. At 3- loops the ar{s) coupling has 
an infrared fixed point whereas at higher orders the behavior of the coupling is not known. 
Estimates of the 4-loop term indicate both an infrared fixed point as well as a divergent 
behavior. 

Recently it has been argued by Howe and Maxwell [49] that the effective charge q;r(s) 
has an infrared fixed point to all orders in perturbation theory. Given the commensurate 
scale relation between ar and aR this should then also be true for ar{s). A simple example 
of the mechanism behind such a perturbative infrared fixed point is given by the effective 
"time-like" one-loop coupling aes{s) which agrees well qualitatively with the empirically 
determined coupling. Another mechanism which also gives an infrared fixed point in qual- 
itative agreement with the data is given by the perturbative gluon condensate dynamics. 
These simple examples show that indeed the empirical behavior ar{s) is consistent with an 
infrared fixed point. 

As we have discussed in the introduction, effective charges defined from different observ- 
ables are related in perturbative QCD to each other at leading twist through commensurate 
scale relations [10]. The perturbative coefficients which appear in these relations are identi- 
cal to those in conformal QCD; i.e. the theory defined with a zero f3 function. The effects of 
the nonzero /3 function of physical QCD are absorbed into the scale of the effective couphngs. 
Examples are the relations between and and the generalized Crewthcr relation, which 
relates the effective charge defined from the Bjorkcn sum rule for gi{x,Q'^) to the Adler 
coupling defined from e+e^ annihilation through a geometric scries [12]. 

Since ar{ml) has effectively constant behavior at small mass scales, the commensurate 
scale relations between ar and other physical couplings imply that at leading twist all 
effective charges have a similar fixed-point behavior in the infrared regime. In particular, 
our results are consistent with the infrared freezing of aR at low s obtained by Mattingly 
and Stevenson [47]. 

The commensurate scale relations [10] can of course also be used at larger scales in order 
to compare experimental results from different observables. For example, the ar effective 
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charge can be evolved to the commensurate scale of another observable which one wants 
to compare with. This way the results can be compared directly, free of theoretical con- 
ventions, instead of by first extracting q;ms(M|) and then comparing. Beyond leading twist 
one can also include power-corrections in the commensurate scale relations and determine 
non-perturbative parameters in a standard way. Ultimately one could also envision to use 
a physical effective charge in the dispersive approach. As we have shown, the four-loop 
evolution of with ^4^^ = 25 gives a good description of the data down to s ~ 1 GcV^ 
and could thus be used as an effective model for the coupling in the low-energy regime. 

The near constancy of the effective QCD coupling at small scales helps explain the 
empirical success of dimensional counting rules for the power law fall-off of form factors and 
fixed angle scaling (for a review see ref. [51]). As shown in ref. [52], one can calculate the 
hard scattering amplitude Th for such processes [53] without scale ambiguity in terms of 
the effective charge or aR using commensurate scale relations. The effective coupling is 
evaluated in the regime where the coupling is approximately constant, in contrast to the 
rapidly varying behavior from powers of predicted by perturbation theory (the universal 
two-loop coupling). For example, the nucleon form factors arc proportional at leading order 
to two powers of evaluated at low scales in addition to two powers of l/g^; The pion 
photoproduction amplitude at fixed angles is proportional at leading order to three powers 
of the QCD coupling. The essential variation from leading-twist counting-rule behavior then 
only arises from the anomalous dimensions of the hadron distribution amplitudes. 

The large magnitude that we find a-r ~ 0.9 ± 0.1 also implies a substantially larger nor- 
malization for the pion form factor and other exclusive observables than estimates based 
on a canonical value as ~ 0.3. The corresponding larger value of a" associated with the 
exchange of n giuons in the hard scattering amplitude could eliminate much of the discrep- 
ancy between data and previous estimates of the PQCD predictions for the normalization 
of exclusive hard scattering amplitudes. 



APPENDIX A: THE /3-FUNCTION IN THE MS AND ol^ SCHEMES 

The perturbative expansion of the {3 function is given by, 

= -Poalis) - ^al{s) - ^,at{s) - ^,al{s) (Al) 

where as{s) — Q;s(s)/(47r). The first two terms in the /3-function [54-58], 

2 

A = 102 - ynf, 

are universal at leading twist whereas the higher order terms are scheme dependent. In the 
MS scheme the first two scheme dependent coefficients are known [59-61]: 

2857 5033 325 2 
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In case of the scheme, 13^,2 is known exactly whereas for (3t,3 there are only estimates. 

Since the three- and four-loop coefficients in the /3-function are known in the MS-scheme, 
the corresponding coefficients (or estimates thereof) can be obtained from the perturbative 
expansion of a^. in the MS-scheme. Starting from the perturbative expansion of the Adler 
D-function for the photon vacuum polarization: 



^, , .2 dn^(-s) 
D(-s) = -^-K^s- ^ ' 



ds 



n=l 



TT" 



with the known coefficients [62-66] 
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one can calculate the following expression for where rzf = 3 is assumed for the 
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where, as is customary in the r literature, we use the following notation for the D„ with 
f = {u,d,s}: 



= L>i = 1, 
K2 = L'2(nf = 3) 

'MS _ T-iMS/ 
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The last coefficient above is only an estimate [30] but it is often used to evaluate the uncer- 
tainty resulting from the missing higher order terms. 

Finally, the resulting non-universal /3-function coefficients in the a^-scheme are then 
given by: 



/26TX3 \ 
+ + 472 (3 + 600 C5 - 1296 C| j l^o 



79813 
16 



+ 6552 Ca + 5400 Cs - 243 tt^ - 11664 Q ^ -789, 



/3..3 = Pf' + + + " ^^^^ + ^^^^ ^0 
/845 5 2\ ^ ^2 f 12673115 ,,,,,, ^ 
-/5o + (y44 - 3^ j /^i /^o + ( 2^ + 555556 Cs 



432 

- 179400 Cs + 129600 Ca Cs + 101952 Q - 186624 C| + 128 Kf^^ (3q 



585179735 

+ 4820288 Ca - 1614600 Cs + 1166400 Ca Cs 



144 _ _ 

+1000512 C|- 8640 TT^- 1679616 C| + 1152 -46776 + 1152 Xf^. (A7) 

Our results agree with the numerical results published in 
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